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Abstract 

We report magnetic field spectroscopy measurements in carbon nanotube quantum dots ex- 
hibiting four-fold shell structure in the energy level spectrum. The magnetic field induces a large 
splitting between the two orbital states of each shell, demonstrating their opposite magnetic mo- 
ment and determining transitions in the spin and orbital configuration of the quantum dot ground 
state. We use inelastic cotunneling spectroscopy to accurately resolve the spin and orbital con- 
tributions to the magnetic moment. A small coupling is found between orbitals with opposite 
magnetic moment leading to anticrossing behavior at zero field. 

PACS numbers: 73.22.-f, 73.22.Dj, 73.23.Hk, 73.63.Fg 
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The remarkable electronic behavior of carbon nanotubes (CNTs) originates from a par- 
ticular combination of the symmetry properties of the graphene band structure and the 
quantization of momentum imposed by periodic boundary conditions along the nanotube 
circumference □ D . This symmetry results in a four-fold degenerate shell structure in the 
energy spectrum of CNT quantum dots (QDs). In early experiments on CNT QDs P, Q], 
however, this symmetry was not observed, presumably due to the presence of defects. Im- 
provements in the quality of CNTs and advances in nanofabrication techniques have enabled 
the recent observation of four-fold degeneracy in the spectrum of CNT QDs 0, (| . An inter- 
esting effect related to the symmetry of the graphene band structure is the modulation of 
the energy gap in CNTs when placed in a parallel magnetic field, B. This effect, predicted 
early on (tJ , has only been recently observed in CNT QDs 0, 0] . These studies, however, 
did not show evidence of four-fold symmetry and the link between the energy spectrum and 
the .B-evolution of the QD states was not established. 

In this Letter, we report -B-dependent electronic transport spectroscopy measurements 
on CNT QDs exhibiting four-fold shell structure. We show that: (i) each shell consists of 
two orbitals with opposite magnetic moment; (ii) the splitting of the orbital states with B 
accounts for all the observed transitions in the spin and orbital configuration of the CNT 
QD; (iii) a weak coupling exists between orbitals with opposite magnetic moment resulting 
in level repulsion at B=0; (iv) Zeeman and orbital contributions to the electron magnetic 
moment can be distinguished by inelastic-cotunneling spectroscopy. 

The electronic structure of CNTs can be derived from the two-dimensional band structure 
of graphene. The continuity of the electron wave function around the CNT circumference 
imposes the quantization of the wave- vector component perpendicular to the CNT axis, k±. 
This leads to a set of one-dimensional subbands in the longitudinal direction Due to 
symmetry, for a given subband at k±=-k a there is a second degenerate subband at k±=k Q . 
Figure la shows in black solid lines the schematic ID band structure of a gapped CNT near 
the energy band gap. Both valence and conduction bands have two degenerate subbands, 
labelled as "+" and "— ". Ajiki and Ando |T| predicted that the orbital degeneracy should 
be lifted by a magnetic field parallel to the CNT axis (Fig. la). This effect can be understood 
by noting that, due to clockwise and anti-clockwise motion around the tube, electrons in 
degenerate "+" and "— " subbands should have opposite orbital magnetic moments, [i or b- In 
the case of finite-length CNTs, a discrete energy spectrum is expected from size quantization. 
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The level spectrum of a CNT QD can then be described as two sets of spin-degenerate 
levels, and with n= 1,2, 3,... (see Fig. la). In the absence of inter-subband mixing, 
E^ = E^ at 5=0, and a four-fold degenerate shell structure is expected for every n. 
(Below we show that a finite coupling can exist, resulting in a small orbital splitting even 
at B=0.) 

The four-fold shell filling emerges in a measurement of the linear conductance, G, versus 
gate voltage, Vq. This is shown in Fig , lb for a QD device fabricated from a metallic 
nanotube with a small band gap P. lid. Ill|. G exhibits Coulomb blockade oscillations 3] 
corresponding to the filling of the "valence" band of the CNT. From left to right, electrons 
are consecutively added to the last three electronic shells, n=3, 2 and 1, respectively. The 
shell structure is apparent from the V^-spacing between the Coulomb oscillations. The 
addition of an electron to a higher shell requires an extra energy cost corresponding to the 
energy spacing between shells. This translates into a larger width of the Coulomb valley 
associated with a full shell Q|. The first group of four Coulomb peaks on the left-hand 
side of Fig. lb („=3) are str„ng. y overlapped due to a large tunnel eoup.ing to the leads 
and Kondo effect |13J. The coupling decreases with Vq, becoming very small near the band 
gap, which lies just beyond the right-hand side of the V^-range shown. Due to this small 
coupling, the Coulomb peaks associated with the last two electrons in n=l are not visible 
(see Fig 2b). 

The shell structure breaks up at finite B (Fig. 2a). In each group of four Coulomb peaks, 
the first (last) two peaks shift towards lower (higher) Vg- This behavior demonstrates the 
strong independence of the orbital levels described in Fig. la. The magnetic field shifts 
the "— " orbital levels down in energy, while the "+" orbitals are shifted up due to their 
opposite /iorfe- Consequently, the addition of the first (last) two electrons to a shell results 
in a pair of Coulomb peaks shifting toward lower (higher) Vq. For each shell, 

l^orb can 

be extracted from the shift, AVcin), in the position of the corresponding Coulomb peaks. 
Neglecting the Zeeman splitting, we use the relation eaAVc(n) = \{i or b(n)costpAB\, where 
AB is the change in B, if is the angle between the nanotube and B, and a is a capacitance 
ratio extracted from non-linear measurements. The values obtained (0.90, 0.80 and 0.88 
meV/T, for n=l, 2 and 3, respectively) are an order of magnitude larger than the electron 
spin magnetic moment (1/2(7/2^=0.058 meV/T f or q =2), and in good agreement with an 
estimate of fi or b based on the nanotube diameter |l4|. 
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The strong B- dependence of the orbital states induces changes in the orbital and spin 
configuration of the QD. These are reflected as kinks in the evolution of the Coulomb peaks 
with B (Fig. 2b). Remarkably, a fully consistent description of the S-dependent energy 
spectrum and the ground state spin/orbital configuration can be obtained through a careful 
analysis of all the kinks in Fig. 2b, as illustrated by the diagrams in Fig. 3. As an example, 
we examine the non-trivial evolution of Coulomb peak CC (notation defined in Fig. 3 
caption). Segment CCi separates the triplet state in region II from the spin 1/2 state in 

(2) 

region III. The incoming electron tunnels into the E_ orbital, with spin down, so the slope 
of the CCi segment is u —^ orb + 1/2^//_b", as noted underneath the corresponding arrow. 
At BxCi, a triplet-singlet transition occurs. Therefore C1C2 separates the singlet state in 
region II from the spin 1/2 state in region III. Now the incoming electron tunnels into the 

(2) (2) 

E y + orbital state with spin up, so the slope of C1C2 is a ^ J b — 1/2^b". Interestingly at 
C2, a kink related to an inter-shell orbital crossing occurs. C2C also separates the singlet 
state in region II from the spin 1/2 state in region III (as C1C2), but the incoming electron 
tunnels into the E_ state and with spin up, so the slope changes direction and has a value 
" — ftorb ~ I/^StW- The rest of the diagrams can be followed in a similar manner. 

Note that kinks in Fig. 2b are connected by conductance ridges crossing the Coulomb 
valleys. The enhancement of G at these ridges is due to Kondo effects of different origins. At 
BxCi, DxEx, and FxGi the Kondo effect arises from singlet-triplet degeneracy 3|. At AB, 
CD, and EF an enhanced Kondo effect is observed in relation to orbital degeneracy The 
Kondo ridges at C2D2 and E2F2 are due to the recovery of orbital degeneracy between E_ 
and E { + ] pjj. Note that, CIS db result of electron-hole symmetry, region III (three electrons in 
shell n—2) and region V (one electron in shell n=l) have a certain degree of mirror symmetry, 
both in terms of the slope of the Coulomb peaks' evolution with B and the Kondo ridges. 

The data shown so far have been explained in terms of a I?-induced splitting of orbital 
degeneracy, as if the two orbital states of every shell were indeed degenerate at B=0. A 
small zero-field orbital splitting may in fact exist and be masked by the Kondo effect at AB, 
CD, and EF. In order to investigate this possibility, we considered a different device with 
a much smaller coupling to the leads and hence much weaker Kondo effect. This device 
also exhibits four-fold periodicity in the Coulomb peaks' pattern. Fig 4c shows a Coulomb 
diamond corresponding to one electron in a shell at 5=80mT 16]. Inside the diamond, 
single electron tunneling is suppressed and transport occurs via higher order cotunneling 
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processes. A sharp increase in the differential conductance, dl/dV, is observed at a bias 
\V\ = Vi n ~ 190/iV, denoting the onset of inelastic cotunneling (IC) 0,Q,Q. The IC 
transition takes place between the two spin-degenerate orbital levels of the same shell thereby 
indicating the existence of a finite splitting at B=0. Before discussing the -B-dependence of 
the IC edges we note that a weak Kondo peak is also present at V=0 (top inset in Fig. 4a). 
This Kondo effect arises from the single-electron occupancy of the spin degenerate orbital 
ground state. 

At finite B, both the Kondo peak and the IC edges split due to Zeeman spin splitting. 
This is shown in Fig. 4a, where dl/dV is plotted vs (V,B) for Vq at the center of the 
Coulomb diamond In order to identify the dl/dV steps more clearly, Fig. 4b shows 
the numerical derivative of the dl/dV plot in Fig. 4a (i.e. d 2 I/dV 2 vs V and B). IC steps 
in Fig. 4a turn into peaks (V>0) or dips (V<0) in Fig. 4b. The zero-bias Kondo peak 
evolves into two dl/dV steps at V = g/jbsB/e (g — 2). These correspond to IC processes in 
which the spin state of the QD is flipped, i.e. from |— , f > (ground state) to |— , |> (excited 
state). Each of the two dl/dV steps associated with inter-orbital IC splits by gusB/e and 
they move further apart due the increasing orbital splitting, 2fi or i,Bcos(p (<p = 33°). We 
estimate 2[i OTb ~ 350/ieV/T i.e. ~ 3 times the Zeeman splitting [20]. The inner inter- 
orbital steps correspond to IC from |— , |> to |+,|>, and evolve with a slope ±2/i or &/e. 
The outer inter-orbital steps correspond to IC from |— , |> to |+,|>, and evolve with a 
slope ±(2/i or 6 + gixsB)/e. The six steps ("Zeeman", "orbital" and "orbital+Zeeman") can 
be seen in the bottom inset to Fig. 4a. Such separation between the orbital and Zeeman 
contributions to the magnetic moment of electrons in CNTs has not been shown before. 
The evolution of the outer IC peaks is non-linear at low B, indicating an anti-crossing 
between the "+" and "— " orbital levels. Such IC spectrum can be readily modelled using a 
Hamiltonian that includes an inter-orbital coupling besides orbital and Zeeman terms. The 
corresponding energy eigenstates are: E = ±^/(5/2) 2 + (/i orb B) 2 ± l/2gn B B (the 4 possible 
sign combinations). The IC spectrum calculated with this simple model (Fig. 4d) clearly 
accounts for the experimental data. The non-linear evolution of the orbital splitting with 
B constitutes direct evidence that the so-called " subband level mismatch" , usually denoted 
by 5 [5] , is due to a small, but finite, quantum mechanical coupling between the two orbital 
subbands in carbon nanotubes. 

We finally comment on the reproducibility of the results shown. G(B, Vq) patterns with 
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a strong orbital contribution to the magnetic moment of electrons (similar to that in Fig. 
2b) have been measured in all devices (five from three different fabrication runs) where 
four- fold shell filling was observed |2jj. Our study demonstrates that the spin and orbital 
configuration of CNT QDs can be understood and controlled by means of a magnetic field. 
By varying B and the angle nanotube-5, researchers have two "semi- independent" knobs 
(B controls the Zeeman splitting and the angle nanotube-I? controls the orbital splitting), 
which will prove very useful in a variety of experiments with CNT QDs, such as the study 
of the Kondo effect in degenerate systems or the interaction between orbital states at high 
B. 
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FIG. 1: (a) Schematic band structure of a CNT near its energy gap. Black lines represent the 
one-dimensional energy dispersion relation, E(ku), at B=0 (fcn is the wave vector along the CNT 
axis). The valence (conduction) band has two degenerate maxima (minima). Size quantization 
in a finite-length CNT results in a set of discrete levels with both spin and orbital degeneracy. 
The degeneracy is lifted by a magnetic field parallel to the CNT. The ID subbands (and the 
corresponding levels) at finite B are represented by red and blue dotted (solid) lines. Only the 
orbital splitting of the energy levels is shown in this figure, (b) Linear conductance, G, vs gate 
voltage, Vg taken at T=8K. Inset: Device scheme. 
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(a) < n = 3 > + n = 2 » <--n = 1--- 




FIG. 2: (a) G vs 5 on a color scale at T=0.34K for the Vq range shown in Fig. lb (red indicates 
high G and dark blue low G). (b) Zoom- in of (a). The green dashed lines highlight the evolution 
of the Coulomb peaks with B. They are labelled as AA', BB'... and FF'. These divide the plot 
in different Coulomb blockade regions indicated by the number of electrons in the last two shells 
(white numbers to VI). The high-G regions (indicated by yellow dashed lines) in between Coulomb 
peaks are due to Kondo effect. Orange numbers indicate the spin in each region. On the right side, 
the G is multiplied by 20, so that the triplet-singlet transition is clearly seen along FiGi. 
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FIG. 3: Diagrams representing the orbital and spin configuration in the different regions of Fig. 
2b. Each row follows the S-evolution of a given Coulomb peak. Each section in a given row shows 
two diagrams (separated by a double arrow). These represent the ground state configuration in 
the two regions separated by the corresponding segment of the Coulomb peak evolution with B. 
In the diagrams, different colors refer to different shells: orange, blue and black for n=l, 2 and 3, 
respectively (in our discussion here n=3 is always filled). The last added electron is displayed in 
red, all the others in black. Note that, for each diagram, the final state is the same as the initial 
state for the diagram immediately below (as indicated by the connected green dashed circles). 
We use solid (dotted) lines to represent levels with positive (negative) fJ, rb, shifting down (up) 
with B (Zeeman splitting is neglected because it is an order of magnitude smaller than the orbital 
splitting). The slope corresponding to the B-evolution of the Coulomb peaks is also indicated under 
the double arrow. Top-right inset: qualitative energy spectrum of the CNT QD as a function of B 
(Zeeman splitting neglected). 
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FIG. 4: (a) Color-plot of the differential conductance, dl/dV , vs bias, V, and B, measured in the 
center of the Coulomb diamond (see (c)) at T=30mK. The yellow dashed lines indicate the traces 
shown in the top and bottom insets. Insets: (top) dl/dV trace taken at 5=80mT, showing the 
onset of inter-orbital IC and a small zero bias peak due to ordinary spin 1/2 Kondo effect. The 
vertical axis scale spans from 0.02 to 0.08 e 2 /h. (bottom) Same as top inset, but at B=0. 7T, 
showing the six IC steps, (b) Numerical derivative of the dl/dV plot in (a). The two inner lines 
result from Zeeman splitting of the Kondo peak at B=0. The outer lines represent the S-evolution 
of the spin-split orbital levels, (c) dl/dV vs V and Vg, for a single electron in a shell at i3=80mT. 
(d) Calculated 5-dependence of the IC spectrum for a single electron in a spin degenerate level 
for two coupled orbitals. Red (blue) lines indicate upwards (downwards) steps in dl/dV with 
increasing V. 
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